February 6, 2008 23:23 



WSPC - Proceedings Trim Size: 9in x 6in como2006 



FLUCTUATIONS AND PATTERN FORMATION IN FLUIDS 
WITH COMPETING INTERACTIONS 

A. IMPERIO 

CNISM, Dipartimento di Fisica, Via Celoria 16, 20133 Milano, Italy 
E-mail: alessandra.imperio@mi.infn.it 

D. PINI and L. REATTO 

Dipartimento di Fisica, Universita di Milano, Via Celoria 16, 20133 Milano, Italy 



One of the most interesting phenomena in the soft-matter realm consists in the 
spontaneous formation of super-molecular structures (microphases) in condi- 
tion of thermodynamic equilibrium. A simple mechanism responsible for this 
self-organization or pattern formation is based on the competition between 
attractive and repulsive forces with different length scales in the microscopic 
potential, typically, a short-range attraction against a longer-range repulsion. 

We analyse this problem by simulations in 2D fluids. We find that, as 
the temperature is lowered, liquid-vapor phase separation is inhibited by the 
competition between attraction and repulsion, and replaced by a transition to 
non-homogeneous phases. The structure of the fluid shows well denned signa- 
tures of the presence of both intra- and inter-cluster correlations. 

Even when the competition between attraction and repulsion is not so 
strong as to cause microphase formation, it still induces large density fluctua- 
tions in a wide region of the temperature-density plane. In this large-fluctuation 
regime, pattern formation can be triggered by a weak external modulating field. 

Keywords: Pattern formation; Competing interactions; Numerical simulation; 
Structure factor; Liquid-vapor equilibrium. 



1. Introduction 

Pattern formation is a fascinating phenomenon widely diffused in nature, 
from the stripes of zebras, to sand-piles, to the arms of a spiral galaxy. If na- 
ture is used to developing patterns, the zoology of morphogenesis observed 
in laboratory experiments is equally wide. Flux of matter and/or energy in 
open systems can be addressed as a source of modulation far from thermo- 
dynamic equilibrium, such as in the Turing patterns 1,2 and the convective 
cells of the Rayleigh-Bcnard instabilities. 3 ' 4 
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Self-organization phenomena, however, can set in also in conditions close 
to thermodynamic equilibrium. From this point of view, soft-matter sys- 
tems are a benchmark, in which to tune the interparticle potential and to 
study pattern formation. Ferrofluids and Langmuir monolayers are typi- 
cal two-dimensional (2D) experimental systems in which the formation of 
micro-domains such as stripes, droplets, bubbles, rings, uniformly spread on 
the surface or, conversely, arranged into highly ordered super-lattices, can 
be found. 5-8 Ferrofluid films are colloidal dispersions of magnetic nanopar- 
ticles in a carrier liquid. In these systems each particle is small enough to 
constitute a single magnetic domain with its own dipole moment. Langmuir 
monolayers are often formed at the air-water interface, providing an ideal 
smooth substrate. They arc commonly formed by amphiphilic molecules, 
which consist of a hydrophilic "head" and a hydrophobic "tail" . If the den- 
sity of the layer is sufficiently high, the tails can orient so that the molecular 
dipoles align and repel each other. 9 ' 10 

Amphiphilic molecules can give rise to microphases also in bulk 3D 
systems, if their concentration is sufficiently high. In this case the forma- 
tion of micelles occurs, which are also called "association colloids" because 
their size is similar to that of the colloidal particles (lnm-1/im). Further 
examples of microphases involve block copolymer systems, 11 ' 12 colloidal 
suspensions 13 ' 14 and liquid crystals. 15 

In colloidal systems, which this work is mainly concerned with, the 
competition between different forces acting on different length scales is often 
invoked to explain the formation of particle micro-domains. Within this 
frame, the effective potential is made up of a short-range attraction, favoring 
the condensation of particles, plus a longer-range repulsion, limiting cluster 
growth. 16 The short-range attraction can be due to depletion or dispersion 
forces, while the longer-range repulsion can be due to dipolar interactions 
or to weakly screened Coulomb forces. The regime in which the repulsion 
is sufficiently strong to induce pattern formation is discussed in Sec. 2. 
The effect of decreasing the strength of the repulsion so as to avoid stable 
pattern formation is considered in Sec. 3. Conclusions and outlooks arc in 



2. Bidimensional pattern formation 

In this section we study a 2D fluid, in which the particles interact via a 
pairwisc spherically symmetric potential of the form: 



Sec. 4. 




(1) 
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with 

U(r) = --^r ex P(-f ) + exp("f ), (2) 

in which r is the inter-particle distance and a the hard disk diameter. The 
subscript letters a and r stand for "attraction" and "repulsion" . The values 
of the potential parameters, for the cases analysed in this section, are: R a = 
lcr, R r = 2a and e a = e r = 1. Hereafter every physical quantity is meant in 
reduced units. Energies and temperatures are in units of Uo = Ujj(r = a), 
lengths in units of a, number densities p in units of a 2 , wave vectors in 
units of er -1 . 
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Fig. 1. Main frame: the pair-potential as a function of the interparticle distance. Inset: 
Fourier transform of the pair-potential. 

The pair potential is characterized by the presence of a repulsive hump 
as shown in Fig. 1. The Fourier transform of the potential is characterized by 
a minimum occurring at a non zero wave vector k m . When the longer-range 
repulsion is absent (e r = 0), such a minimum occurs at k m = 0. Roughly 
speaking the position of the repulsive hump controls the typical cluster size, 
while the wave vector at which the minimum of the Fourier transform of 
Uff(r) occurs governs the typical inter-cluster distance or pattern period 
P ~ 2rr/fc m . 17 

This model has been studied mainly via numerical simulations in the 
canonical ensemble (NVT, N number of particles, V volume, T tempera- 
ture). We adopt the parallel tempering (PT) 18-20 technique, which is es- 
sentially a Monte Carlo technique, but applied to an extended ensemble, 
involving many different temperatures. 

At low temperatures, we observe the formation of micro-domains whose 
shape depends on the mean density of the system. In particular, at low 



e=1 £ a -1 R=2 R a =1 
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density, we have the formation of small circular droplets. On increasing the 
density, the droplets enlarge and order onto a triangular super-lattice. For 
densities larger than ~ 0.3, the domains morphology switches from circular 
to striped. For densities greater than ~ 0.55, we observe the formation of 
bubbles of gas spread into a liquid phase; this case is somehow specular 
to that of the super-lattice of liquid droplets. The different pattern mor- 
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Fig. 2. Static structure factor computed along the y direction at p = 0.4 and T = 0.5. 
Two series of peaks can be identified, which are related respectively to the correlations 
between different micro-domains (k < 3) and within the same micro-domain (fc > 3). 



phologics affect the features of the static structure factor S'(k), which is 
calculated by explicit evaluation of the expression: 

S(k) = N- 1 < PkP _ k >= 

N N 

= N ^ < E cos ( k • r *)) 2 + Gl sin ( k • r *)) 2 >> ( 3 ) 

i i 

where <> indicates the ensemble average. In the structure factor profiles, 
two different series of peaks (Fig. 2), corresponding to different length scales 
in direct space, can be identified: 21 the peaks at short wave length (k < 3) 
are related to the spatial modulation of the density in the microphases, 
that is they represent the cluster-cluster correlations and we will refer to 
them as "modulation peaks"; the peaks at higher wave vectors (fc~3) are 
connected to the particle-particle correlations within clusters. 

In Fig. 3 we have plotted both the snapshots of the systems at different 
densities (left panels) and the corresponding S(k) at short wave vectors. 
For low densities, such as p = 0.05, the maximum values of S(k) are rather 
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large (~ 18) but the peak is independent of the direction of k. This spherical 
symmetry of S(k) suggests that the droplets are disordered, so that the 
whole structure resembles that of a liquid, in which the droplets are the 
"effective" particles. When the micro-domains order onto a super-lattice, 
the static structure factor exhibits high Bragg peaks, whose symmetry is 
sixfold for the triangular lattice of droplets or bubbles, while it is twofold for 
the striped phase. On increasing the temperature, 5(k) becomes isotropic 
at all densities, showing a ring of small amplitude at k m , untill it totally 
disappears when the system becomes homogeneous. 

3. The weak long-range repulsion regime 

If the strength e r of the repulsive long-range potential in Eq. (2) is suffi- 
ciently small, the system displays a standard liquid-vapor phase transition 
with a critical temperature, which is a decreasing function of e r . The state at 
which a fluid, subject to competing interactions, stops to undergo a liquid- 
vapor transition favoring a modulated phase is called Lifshitz point (LP). 22 
We have used liquid-state theory in 3D to investigate how the liquid- vapor 
transition is affected by the proximity to a LP. The model potential which 
we used is similar to that of Eq. (1), but we have replaced the simple ex- 
ponentials with Yukawas for reasons of computational ease: 

w--^«p(-£) + i *«p(-£). w 

This potential has been studied by the self-consistent Ornstein-Zernike ap- 
proximation (SCOZA). 23 SCOZA has proven to be an accurate and robust 
approach, which lends itself to a semi-analytical formulation for potentials 
consisting of a superposition of Yukawa tails, like that considered here. 
The ranges of the attractive and repulsive contributions in Eq. (4) were 
set at R a = la, Rb = 2a as before, while the repulsion-to-attraction ra- 
tio A = e r /e a has been varied so as to get as close as possible to the LP. 
Since SCOZA is not able to deal with the formation of non-homogenoeus 
structures such as those described in the previous section, this ratio cannot 
be increased at one's will: above a certain value of A, the theory breaks 
down before a liquid-vapor phase transition occurs. We take this as the 
limit beyond which liquid-vapor phase separation ceases to be stable, and 
is replaced by microphase formation. Simulations of fluids with competing 
interactions have indeed shown that micro-domains occur also in 3D. 24 

As the limit of stability of the liquid-vapor transition is approached, 
there are some interesting effects that come along. Figure 4 compares the 
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Fig. 3. Left panels: snapshots of the system for different mean densities. Right panels: 
static structure 5(k) = S(k x ,k y ) computed at short wave vectors k. Top to bottom: 
p = 0.05 T = 0.3, p = 0.20 T = 0.5, p = 0.40 T = 0.5, p = 0.60 T = 0.5. 

liquid-vapor equilibrium curves for a fluid with purely attractive tail po- 
tential (corresponding to A = 0) and one close to the LP (corresponding to 
A = 0.059). When competition between attraction and repulsion is present, 
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Fig. 4. Liquid-vapor equilibrium curve of the fluid with interparticlc potential U(r) 
given by Eq. (4) for A = € r /e a = (dashed line) and A = 0.059 (solid line). Density p 
and temperature T have been rescaled by their critical values p c , T c . The filled region 
is the supercritical locus for A = such that the reduced compressibility Xrod is larger 
than 10 2 . The hatched region is the corresponding locus for A = 0.059. 

the equilibrium curve is much natter, and resembles that expected for a 
first-order transition. Moreover, the supercritical region where the isother- 
mal compressibility of the fluid is large is greatly enhanced. This is shown 
again in Fig. 4, where the locus of supercritical states such that the com- 
pressibility exceeds one hundred of its ideal-gas value is shown, bith with 
and without competition. Since the compressibility of a fluid can be taken 
as a measure of its density fluctuations, it appears that, when competi- 
tion is present, the fluid manages to achieve large density fluctuations in 
a substantial domain of thermodynamic states without resorting to phase 
separation. Further insight into this scenario is gained by inspection of the 
correlations. In Fig. 5 we have plotted the radial distribution function g(r), 
whose Fourier transform is straightforwardly related to the structure factor 
mentioned in Sec. 2, again for the cases A — and A = 0.059. In both cases, 
the thermodynamic states considered are close to the liquid-vapor critical 
point, and correspond to similarly large values of the compressibility. How- 
ever, the strength of the correlations appears to be very different in the two 
cases. At first this may seem somewhat puzzling, since the compressibility 
of the fluid is related to the integral of g(r) by the compressibility rule 



Xrcd (the "reduced" compressibility) being the ratio between the compress- 
ibility and its ideal-gas value. However, upon closer inspection, one finds 
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Fig. 5. Radial distribution function of the fluid with intcrparticlc potential U(r) given 
by Eq. (4). Dashed line: A = 0, p = 0.275, T = 2.496. Solid line: A = 0.059, p = 
0.246, T = 1.920. Densities have been set to their critical values, and temperatures 
have been adjusted so as to give similarly large values of the reduced compressibility, 
Xrcd = 3.8 X 10 5 and x rc d = 2.2 X 10 5 . The inset shows the large-r limit of g(r). 



that Eq. (5) is consistent with the behavior shown in Fig. 5 since (see inset) 
the g(r) of the fluid with purely attractive interactions dominates for large 
enough r. In other words, when faced with the task of achieving large den- 
sity fluctuations, the fluid with competing interactions "chooses" to enhance 
particle-particle correlations at short and intermediate distances compared 
to the situation found in the absence of competition, while depressing them 
at very large distance. Such a strong enhancement of correlations at inter- 
mediate lengthscales indicates that the density fluctuations responsible for 
the large values of Xrcd must be traced back to the occurrence of delocalized 
clusters of strongly correlated particles. These can be considered as the pre- 
cursors of the non- homogeneous phases discussed in Sec. 2, that eventually 
set in at higher values of A. This conclusion is further supported by an 
analytical study of the decay of the correlations, which shows that, close to 
a LP, there are two characteristic lengths at play: one of them reduces to 
the usual correlation length, and hence diverges as the liquid-vapor criti- 
cal point is approached, while the other one remains finite, but still much 
larger than the particle diameter, and is related to the size of the region of 
enhanced correlations. 

The behavior described above also suggests that in this region the sys- 
tem should be extremely sensitive to the presence of an external modula- 
tion. We have studied this condition in 2D via simulations, adopting the 
interaction potential of Eqs. (1), (2), when the system its close to the LP. In 
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this case the potential parameters (e a = 0.679, e r = 0.22, R a = 1, R r = 2) 
are chosen in such a way that the minimum of the Fourier transform of the 
pair potential is at k m = 0. Therefore, even though the fluid is subject to 
competing interactions, it is no longer characterized by an intrinsic peri- 
odicity. In particular, we have analysed the response of the system to the 
presence of an external agent, such as a potential modulated along the x 
direction. Such a potential induces the formation of parallel stripes, whose 
period is the same as that of the modulating potential. Comparing the re- 
sults with those obtained from a standard Lennard- Jones fluid, subject to 
the same external modulation, we find that our system is much more af- 
fected by the action of the field. In other words, next to the Lifshitz point, 
the fluid experiences large density fluctuations which can be easily driven 
into a stable microphase by the presence of a very weak external potential. 21 

4. Conclusions and outlooks 

In this work we have analyzed the behavior of a fluid system at ther- 
modynamic equilibrium, in which the particles are subject to competing 
interactions. Such a model represents a very simple and effective micro- 
scopic description of the mechanisms which can induce pattern formation 
or self-organization, such as those observed in many experimental situa- 
tions involving colloidal particles both in 2D and 3D. Even though the 
interaction potential is characterized by spherical symmetry, we have ob- 
served the formation of very different micro-domain morphologies such as 
droplets, stripes and bubbles. 

When the long-range repulsion is not sufficiently strong to stabilize 
any pattern formation, there exists all the same a wide region in the 
temperature-density plane, in which the fluid is subject to large density 
fluctuations. In such a region the fluid is very sensitive to an external mod- 
ulating field. 

In soft matter, nowadays, a great effort is devoted to understand not 
only the general mechanisms supporting pattern formation in bulk systems, 
but also the effect of a geometrical confinement due to the presence of con- 
fining walls or external fields. A study 26 carried out in the weak long-range 
repulsion regime of Sec. 3 shows that, because of the tendency towards clus- 
tering caused by the competition, density profiles at wall-fluid and liquid- 
vapor interfaces are considerably more structured than those found in fluids 
with purely attractive interactions. Moreover, both experiments and simu- 
lations show that confinement is a powerful tool to modify and even induce 
new micro-domains morphologies, which in the bulk might not appear at 
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all. This is particularly expected when the extent of the spatial region is 
not commensurate with the intrinsic periodicity of the bulk system. 
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